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ABSTRACT 

We study the effects of gravity waves, or g-modes, on hot extrasolar planets. These planets are expected to 
possess stably-stratified atmospheres, which support gravity waves. In this paper, we review the derivation of 
the equation that governs the linear dynamics of gravity waves and describe its application to a hot extraso- 
lar planet, using HD 209458 b as a generic example. We find that gravity waves can exhibit a wide range of 
behaviors, even for a single atmospheric profile. The waves can significantly accelerate or decelerate the back- 
ground mean flow, depending on the difference between the wave phase and mean flow speeds. In addition, the 
waves can provide significant heating (~ 10 2 to ~ 10 3 K per planetary rotation), especially to the region of the 
atmosphere above about 10 scale heights from the excitation region. Furthermore, by propagating horizontally, 
gravity waves provide a mechanism for transporting momentum and heat from the dayside of a tidally locked 
planet to its nightside. We discuss work that needs to be undertaken to incorporate these effects in current 
atmosphere models of extrasolar planets. 

Subject headings: hydrodynamics — planets and satellites: general — waves — turbulence 



1. INTRODUCTION 

A stably-stratified atmosphere, characterized by a positive 
vertical entropy gradient, can support gravity waves, or g- 
modes. Gravity waves are oscillations which arise from the 
buoyancy of parcels in the fluid. Such waves are readily ex- 
cited by flow over thermal and surface topography, convective 
and shear instabilities, and flow adjustment processes. They 
propagate through the atmosphere both horizontally and ver- 
tically. Gravity waves i n the terrestrial atmosphe r e and ocea n 
are much studied (e.g.. iGossard & Hookdfl97l lGiIilll982l) . 
They have also bee n observed on oth er Solar Sy stem bod- 
ies, su ch as Jupiter (I Young et al.ll 19971) and Venus (lApt et al.l 
[l98(ih . 

In the terrestrial atmosphere, a typical gravity wave has an 
energy flux of approximately 10~ 3 to 10 _1 W m" 2 . Despite 
being small, compared to the total amount of absorbed solar 
flux (^237 W m~ 2 ), gravity waves are responsible for sig- 
nificantly modifying — even dictating — large-scale flow and 
temperature structures. Several well known examples of this 
are the Quasi-Biennial Oscillation, reversal of mean merid- 
ional temperature gradient in the upper middle atmosp here, 
and generation of turbulence (e.g.. lAndrews et aT1ll987h . We 
expect similar effects to be present in the atmospheres (and 
oceans) of extrasolar planets. Moreover, due to the greater 
irradiation and scale heights on them, the acceleration and 
heating effects of gravity waves can be much stronger on hot 
extrasolar planets. 

There has been much interest in modeling the at- 
mospheric circu la tion of extrasola r pla ne ts (e.g., 
Joshi et al.1 [19971 IShowman & Guillol l200l ICho et al.l 
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circulation is crucial for interpreting observations of extraso- 
lar planets, as well as for improving theoretical understanding 
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in general. For this, the role of eddies and waves in transfer - 
ring momentum and heat needs to be addressed dCholl2008l) . 
This has long been recognized in Solar System planet studies 
(e.g.. lLindzedll990tlFritts & Alexanderl2003l) . 

The plan of the paper is as follows. In §|2]we derive the gov- 
erning equation appropriate for linear monochromatic grav- 
ity waves on hot extrasolar planets. We also discuss a sim- 
ple parameterization of the key non-linear process, saturation. 
In addition, we present solutions to the equation for simple 
isothermal atmospheres, with and without shear in the back- 
ground mean flow. In ^3]we extend the calculation to a phys- 
ically more realistic situation, by using background flow and 
temperature profiles derived from a three-dimensional (3-D) 
hot-Jupiter atmospheric circulation simulation. This is the 
first such calculation to have been performed for extrasolar 
planets. Through this, the significant effects of gravity waves 
on hot extrasolar planet atmospheric mean flows are demon- 
strated. In this section, we also discuss a way in which gravity 
waves can transport momentum and heat horizontally — e.g., 
from the dayside to nightside on tidally locked planets. In §|4] 
we discuss the implications of our work for current extrasolar 
planet atmospheric modeling work. We conclude in ^5] 



2. LINEAR THEORY 

2.1. Taylor-Goldstein Equation 

The dynamics of a linear gravity wave is described by the 
Taylor-Goldstein equation (TGE). This eq uation is derive d 
from the full, 3-D hydrodynamics equations (Bat chelorll967l) . 
In this work, we restrict the description to two dimensions and 
neglect rotation: 
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where D/Df = d/dt + u ■ V; u = (w, w) is the flow in the hori- 
zontal and vertical directions (x,z), respectively; p is the den- 
sity; V = (d/dx,d/dz)', p is the pressure; g = (0,-g) is the 
gravity; 9 is the potential temperature; c p is the specific heat 
at constant pressure; T is the sensible temperature; and Q is 
the net diabatic heating rate. Equations (HJ are supplemented 
with the ideal gas law: 



p = pRT, 

where R is the specific gas constant. Note that 



T ,Ei 



(2) 



(3) 



where p R is some reference pressure (here taken to be the pres- 
sure at the lower boundary of the model) and k = R/c p ; 9 is 
related to the entropy s by ds = c p din 9. 

The neglect of the third dimension and rotation requires 
that we restrict our analysis to waves with horizontal scale 
L < U/Cl, where U is the characteristic mean flow speed and 

is the planetary rotation rate. This scale is adequate for all 
gravity waves, except for large-scale tides ( which has been re- 
cently considered by iGu & Ogilvid (120091) V As an example, 
U/Q ~ 10 7 m for HP 2 09458 b, based on U in ho t extraso- 
lar planet simulations of iThrastarson & Choi (120091) . The re- 
sulting value is approximately 1/10 of the planet's radius. In 
addition, g, R, and c p are taken to be constant and Q is spec- 
ified. These restrictions do not mitigate the basic application 
and implications presented in this work. However, for broader 
applications, relaxation of these and other restrictions will be 
considered in future work. 

The variables in equations ([TJ are all expanded as a small 
perturbation about a mean value, which is a function of height 
only: 



£(x,z,t) = Co(z) + Ci(x,z,t). 



(4) 



For the thermodynamic variables, we require that Ci /Co ^ 1 ; 
however, this is not required for the flow variables, u and 
w. The mean state is assumed to be in hydrostatic balance, 
dpo/dz = -pog, and contains only horizontal flow so that 
wo = 0. We also assum e the anelastic approximation (e.g., 
lOgura & Phillipslll962l) . V • (p u) = 0. This implies the fol- 
lowing: 
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where N(z) = [g(dln# /dz)] 1//2 is the Brunt-Vaisala fre- 
quency; H p (z) = poCdpo/dz)" 1 and H p (z) = |p (d/7 /dz) _1 1 
are the density and pressure scale heights, respectively; 
7 = c p jc v is the ratio of specific heats, with c v the specific 
heat at constant volume; c s = (jRT) 1 ' 2 is the speed of sound; 
and, D is the vertical scale of the motion. 



With (@]i and ||5}, we obtain 
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where <&i = p\/po and Oi = 9\/9{). Since the coefficients in 
equation (O are independent of x and t , we assume perturba- 
tions of the form, 



Ci(jc,z,0 = C,{z)exp{i(kx-ujt)}, 



(7) 



where it is understood that the real part is to be taken. This 
leads to the polarization equations: 
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where F = Q/(c p Tq) is the forcing, c = tu/k is the constant 
(possibly complex, see §12. 3b horizontal phase speed, and (c— 
uq) is the intrinsic phase speed. Now, transforming to a new 
variable so that the effect of decreasing density with height is 
compensated, 



w(z) = wexp{-x(z)} 



X(z) = 



2HJQ 



where Zb is z at the bottom, we obtain the TGE: 
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Here, m = m(z) is the index of refraction, and corresponds to 
the local vertical wavenumber. It is given by 
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In equation (fTTT i we have used H p = RTo(z)/g and the prime 
indicates differentiation with respect to z. 

The vertical structure of the perturbations, oscillating in z, 
signify that we are dealing with internal waves. In equa- 
tion ( fTTI ). the key terms contributing to m 2 are the first ("buoy- 
ancy") term and the last ("non-hydrostatic") term. Although 
the other three terms contribute, generally the buoyancy and 
non-hydrostatic terms control whether the wave propagates 
since the flow shear and curvature are small and the scale 
height is large in practice. For large wavelength waves, the 
waves are hydrostatic and the non-hydrostatic term is small; 
then, the buoyancy term dominates. In these cases, as long 
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FIG. 1. — A gravity wave propagating in an isothermal (7b = 1350 K) and 
constant background flow (no = 350 m s _I ) atmosphere. The phase speed of 
the wave c is 100 m s _1 and the horizontal wavelength, 2n/k, is 2500 km. 
The vertical perturbation velocity w (■ ■ ■), horizontal perturbation velocity 

u ( — ), and wave stress r ( ) are shown. The latter is the vertical transport 

of the horizontal momentum. Wave amplitudes, u and w, grow exponentially 
with height, but the wave stress is constant with height, since there is no 
dissipation. The jump in r near z/H p = 1 is caused by the forcing. 

as the atmosphere is stratified (i.e., N 2 > 0) and c uq, the 
wave will propagate vertically. However, for shorter, non- 
hydrostatic waves, it is possible that k 2 > [N/{c-uq)] 2 . In 
these cases, m is imaginary and the wave do es no t propagate. 
This situation is discussed in more detail in 33.21 

Equation dTOb can be thought of as a driven harmonic os- 
cillator equation. When m is real and constant, its solution 
is a simple sinusoid. Since the transformation, equation (0, 
compensates for the exponential decay of density with height, 
w grows (decays) rapidly with height (depth). Here, we have 
dropped the tilde overscript. From hereon we drop all tilde 
overscripts and "1" subscripts for notational clarity. The 
growth can be clearly seen in Figure Q] along with the corre- 
sponding constant stress (vertical transport of h orizontal mo- 
mentum) when there is no dissipation (see 32.31 l. 

If the temperature or the flow varies with height, m is a func- 
tion of height z. If m va ries slowly, the WKB approximation 
dBender & Ors zag 1999) can then be used to obtain 
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Here, A = w(zh)[m{zh)] 1 ^ 2 - As in the constant m case, the 
vertical perturbation velocity is wave-like with upwardly and 
downwardly propagating components; the amplitude of the 



upward component grows with height and the downward 
component decays with depth. However, when m varies 
rapidly, the solution must be obtained numerically. 

At the boundaries we use the radiation condition, selecting 
the upwardly propagating solution at the top boundary z, and 
the downwardly propagating solution at the lower boundary 
Zb- This is achieved using the condition, 
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Here, s = ±1, depending on the signs of the horizontal and 
intrinsic phase speeds and on whether the condition is at the 
upper or lower boundary. Note that condition ( fT3l ) requires the 
WKB approximation to be valid at the boundaries. For exam- 
ple, the boundaries cannot be critical layers, regions where 
(c- Mo) — > 0, since m — > oo approaching such a layer and the 
WKB approximation ceases to be valid. Critical layers will 
be discussed further below. 

As alluded to above, we use Gaussian elimination 
(ICanuto et al.l 12006b in this work to numerically solve the 
TGE. In the cases presented here, 3000 equally spaced levels 
are employed to solve for w. We have checked that this res- 
olution is adequate by performing calculations with 10,000 
levels and verifying convergence at the higher resolution. Ex- 
tensive validations against known analytic solutions, where 
they exist, have also been performed. 

2.2. Polarization Relations and Fluxes 

The solution to the TGE is a wave in the vertical velocity 
perturbation. This can be related to the horizontal and tem- 
perature (sensible and potential) perturbations. Understand- 
ing these are essential for parameterizing the saturation pro- 
cess in the full non-linear situation and in general circulation 
models. The perturbati on quantities can be obtained from the 
polarization equations dHineslll960l) : 
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where we have introduced the geopotential perturbation func- 
tion $ (= p/po)- Equation (fl4h) gives the relationship be- 
tween the vertical and horizontal velocity perturbations. Note 
that the amplitude of u is larger than the amplitude of w, as 
it is scaled by kH p . Note also that the geopotential (pres- 
sure) perturbation varies with the background flow via the 
dependence on the intrinsic phase speed, according to equa- 
tion ([14b). Equations (fl4b) and (fl4h) shows that the potential 
and sensible temperature perturbations, respectively, are both 
7r/2 out of phase with the vertical velocity perturbation. But, 
the phase between u and w varies locally through their depen- 
dence on the background. 

Gravity waves are an efficient means of transporting both 
momentum and energy. The (perturbation) momentum and 
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energy fluxes are simply obtained from equations ( TT4T >: 



T = PqUW 

F x = po$u 
F z = po &w . 



(15a) 
(15b) 
(15c) 



In equations ([TBl l. r is the vertical flux of horizontal momen- 
tum (or, the wave stress); F x and F z are the horizontal and ver- 
tical fluxes of energy, respectively; and, the overline indicates 
an average over a wavelength, 



a/3 = ^Re{a/T}, 



(16) 



where a and /3 are arbitrary complex functions and "*" de- 
notes the complex conjugate. Note that the energy fluxes de- 
pend on the background flow through their dependence on 
po and $. However, as can be seen in Figures Q] and [2j 
the wave stress remains constant, away from the forcing and 
damping regions — e.g., saturation regions and critical layers. 
This is in accordance w ith the second Eliassen-Palm theorem 
(El iassen & Palmlll960l) . which expresses non-interaction of 
the disturbance in the absence of dissipation and forcing. 

In Figure|2j it is important to note that, where the fluxes are 
changing, the wave is interacting with the background flow 
in those regions. This should be contrasted with the behav- 
ior illustrated in Figure [T] where the stress is not changing. 
Changes in the wave stress cause accelerations to the mean 
flow. Correspondingly, changes in the energy fluxes cause 
the temperature of the region to change. The rates of these 
changes are given by 
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(17a) 
(17b) 



In this case, the wave causes the background flow to acceler- 
ate from 350 m s -1 towards 500 m s" 1 . 

In the remainder of the paper, when speaking of vertically 
propagating waves, we consider only waves that propagate 
upwards from the region of excitation. However, it must be 
remembered that downward propagating waves are also gen- 
erated. On a giant planet without a solid surface, those waves 
may not be reflected or absorbed. They can continue to pene- 
trate downward until they encounter a critical level or a con- 
vective region. Or, they are simply dissipated since the am- 
plitudes of the downwardly propagating waves decrease ex- 
ponentially, as already discussed (and as also can be seen in 
Figure |2). The energy and m omentum fluxes are link ed by the 
first Eliassen-Palm theorem (Eliassen & Palm I960): 



F- = (c-ho)t, 



(18) 



which can be derived from equations ( TBI and 1151 . For down- 
wardly propagating waves, we have F z < 0; and, therefore, via 
equations (IT8b and (fTTh). we see that for these waves the de- 
position of momentum also leads to acceleration of the flow 
toward the phase speed of the wave. Downwardly propagating 
planetary scal e gravity waves (i.e., t hermally excited tides) are 
considered bv lGu & Ogilvlel (120091) . 

2.3. Saturation and Critical Layers 

The present work concerns inviscid, linear, monochromatic 
waves. Such waves are infinite in extent and, in principle, can 
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FIG. 2. — As in Figure[T] but with negative vertical wind shear: uq varies 
from 350 m s _1 near the bottom to 650 m s near the top, with linear growth 
in between. Here, c = 500 m s ; hence, (c-uq) — > at z/H„ f» 8, where 
the wave encounters a critical layer and dissipates. In the region just below 
the critical layer, the wave saturates and transmits momentum into the back- 
ground flow, as indicated by the drop in r with height there. 

grow without limit when they propagate upward. This is ob- 
viously not physical. In reality, such waves become unstable 
and saturate. The saturation process can be treated by intro- 
ducing a correction to the solution according to the physical 
criterion, 

f- (Oo + O) < 0, (19) 

dz 

ensuring that a wave does not become convectively unstable 
and remains at neutral stability. When the saturation region is 
identified, 9 is adjusted so that neutral stability is maintained. 
This new value for 9 is then used in equation (fT4b ) to obtain 
the corresponding w in the saturated wave. This value of w is 
then used in equations ( fT4k . b, c), ( fT~5T > and (fTTT i. This satura- 
tion condition acquires the simple form, 

\u\ = c — wo , (20) 



in situ ations where the WKB approximation is valid (iFrittsI 
[19841) . 

As can be seen from the definitions of N, H p , and 9, the 
Brunt-Vaisala frequency can be written as 



dT 



"I 1/2 



(21) 



Since g and c p are essentially constants in the modeled height 
range, depends only on the temperature profile T(z). For 
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isothermal regions, N is a constant. In general, the frac- 
tional change of T with height is small compared to g/(Tc p ) 
throughout the modeled region. Hence, N is nearly constant in 
the entire domain, with a value that is roughly 2.4 x 10~ 3 s -1 . 
The Brunt-Vaisala frequency for our atmospheric profile is 
shown in Figure [3] Furthe r discu ssion of the effects of varia- 
tion in N are presented in §13.1.21 

If the background flow contains shear, it is possible for 
the wave to encounter a critical layer at some height. At 
the critical level, where c = uq, the TGE becomes singular. 
However, t he equation can be solv ed using the method of 
Frobenius (B ender & Orszagl 119991) . from which it is seen 
that the wave is drastically attenuated by the critical layer 
dBooker & Brethertonll 19671) . The amount of attenuation de- 
pends on the Richardson number Ri of the flow, 

N 2 

Ri = j • ( 22 > 

(dwo/dz) 

Figure |2] illustrates a critical layer encountered by a wave. 
If the wave stress has magnitude r below the critical layer, 

it then has magnitude, r • exp j- 27r [Ri-(1 /4)] 1//2 |, after th e 

encounter with the critical layer dBooker & Brethertonll 967b . 
This is a substantial amount of attenuation. For example, in 
our model atmosphere of Figure|2] Ri > 900. Hence, the wave 
is essentially completely dissipated at the critical level, with 
the wave stress falling to practically zero and the momentum 
deposited into the mean flow. Note that, during its approach 
to the critical layer, the wave actually saturates and deposits 
momentum and energy over a finite layer (cf., Figure |2). It is 
important to note that while a wave saturates when approach- 
ing a critical layer, the presence of a critical layer is not re- 
quired for saturation. Saturation is a general process referring 
to wave dissipation by many different mechanisms — e.g., ra- 
diation, conduction, breaking, and turbulence. 

Numerically, when a critical layer is present, we lift the 
phase speed from the real axis by adding a small imaginary 
component: c = c, +/c,, where |c,-/c r < 10~ 3 . This intro- 
duces a small amount of linear damping and ensures that 
the neglected nonlinear terms do not dominate in the re- 
gions where waves become steep and eventually break. Of 
course, adding damping causes the wave-stress to decrease 
with height and the second Eliassen-Palm theorem to be no 
longer valid. However, the effect is small. This can be seen in 
Figure |2 where the wave stress falls negligibly over the layer 
from z/H p w 1 to z/H p w 5. 

3. EXTRASOLAR PLANET APPLICATION 
3.1. Setup 

Our aim in the present paper is to demonstrate several prop- 
erties of gravity waves likely to be important for hot extraso- 
lar planets. For this, we choose HD 209458 b as a paradigm 
planet. This planet has been the focus of many theoretical and 
observational studies, and it is expected to be generic with 
respect to the properties discussed here. The physical param- 
eters that characterize the planet's atmosphere are given in 
Table UJ 

3.1.1. Forcing 

In a stratified atmosphere, gravity waves are readily gener- 
ated by many mechanisms — thermal and mechanical — such 
as the absorption of stellar radiation, convective release of 



TABLE 1 
Parameters Used for HD 209458 b 



Specific Gas Constant 


K 


3523 J kg" 1 KT 1 


Specific Heat at Constant Pressure 




12300 J kg -1 K" 1 


Acceleration Due to Gravity 


g 


10 m s" 2 


Rotation Rate 


n 


2.08xl0~ 5 s _1 



latent heat, storms, flows over topography and coherent lo- 
calized "heat islands", and detonation by impacts. In this 
work, we consider small- and meso-scale thermally excited 
waves, rather than large-scale waves, as already discussed. 
The horizontal wavelengths used in this work are 2500 km or 
less. This is a reasonable range since it is well within the ob- 
served range of gra vity waves on Jupiter dYoung et al.ll 1997k 
iHickev et al.l 120001) . Although not unimportant, we do not 
dwell on the precise nature of the source of the excitations. 
The main interest here is in the propagation and deposition of 
momentum and energy. 

The forcing in equation ( fTOb is simply represented as a 
Gaussian, modified so that it is zero beyond two half-widths 
above and below the center. The center is located at z/H p = 1 
above the bottom of the domain. The half-width is 75 km, 
or ~ 0A5H p . The forcing location and width are chosen so 
that the vertical scale of the forcing is less than the vertical 
wavelength of the waves we present here. We have exten- 
sively explored the effects of varying the location, width, and 
strength of forcing and present this case to illustrate several 
important points. Not surprisingly, the dynamics do depend 
on the chosen parameter values, but the dependence is broadly 
predictable. For example, if the forcing scale is much larger 
than the vertical wavelength, then only a very small amplitude 
wave is emitted from the forcing region, due to cancellations. 

The heating rate, Q/c p , is set to 10~ 3 K s . Note that this 
is a modest value. The forcing corresponds to roughly 300 K 
per rotation of the planet. This is compared to ~100 K per 
rotation at the cho s en loc ation in the circulation model of 
iThrastarson & Choi (120091) . A forcing of ~1 100 K per rota- 
t ion, fo r a similar latitude-longitude locat ion on the planet (see 
33.1.21 ). is used in lShowman et al.l d2008l) . The latter value im- 
plies that, in the absence of motion, the location on the planet 
will cool completely in one rotation of the planet. We stress 
that locally — i.e., scales far below the grid scale of the cur- 
rent circulation models — the forcing could actually be much 
stronger. The actual value is presently uncertain and likely to 
be spatially and temporally variable over the planet. To pro- 
vide a context, for the Earth the heating rate is ~ 2 K per day 
(1 day = 0.29 rotation of HD 209458 b) over large areas; but, 
locally, at the tops of low clouds on the Earth the rate can be 
up to ~ 50 K per day (IWallace & Hobbsll2006l) . 

3.1.2. Background Structure 

Figure [3] shows the mean flow and temperature profiles 
used to obtain much of the results presented in this sec- 
tion. The lower part of both profiles — approximately the 
lower six scale heights — is taken from global circulation sim- 
ulations of the hot extras olar giant planet HD 209458 b by 
IThrastarson & Choi d2009l). using the NCAR Community At- 
mosphere Model dCollins et al.l 120041) . The profile is from a 
point near the equator, slightly away from the substellar point 
(70°E, 10°N). This point was chosen as it is within the equa- 
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FIG. 3. — Sample atmospheric mean flow iiq ( — ), temperature To ( ), and 

Brunt-Vaisala frequency, N (• • • ) profiles of a typical hot extrasolar planet 
HD 209458 b, used in this work. The profile is representative of a region at 
approximately 70° E, 10°N. The planet is a close-in giant planet. The pro- 
files are obtained from a 3-D, global circulation model, up to ~ 10~ 3 bar 
level. Above that level the pro files are simply extended, following loosely 
the observed profiles of Jupiter I Young et al. 2005; Flasar et al. 2004). 

torial jet and the Coriolis parameter / = 2fi sin <fi is not large. 
The temperature profile generally increases with height over 
the lowest 4 scale heights and then becomes isothermal. For- 
tuitously, this provides an opportunity for a loose validation of 
our model: it is very similar to the temperature structure ob- 
served by the Galileo probe in the same region of Jupiter's at- 
mosphere dYoung et a l. 2005). We extend the profile by keep- 
ing the temperature isothermal through the planet's strato- 
sphere and having the thermosphere (beginning of the tem- 
perature inversion near the top) start between 12 and 13 scale 
heights at po « 4 x 10~ 6 bar. This profile has a Richardson 
number of at least 3.4, giving an attenuation of 1.4 x 10~ 5 . 
Hence, any critical layer can be considered to fully dissipate 
the wave. 

The chosen flow profile has two local flow velocity max- 
ima. The upper maximum is extended into a jet with a peak 
at z/Hp « 6. This is similar to the structure of th e jet in 
Jupiter's stratosphere proposed by Flasar et al] ( 120041) . which 
has a peak velocity between the 10~ 2 and 10~ 3 bar levels. We 
then extend the profile further upwards without shear. Note, 
this is in keeping wit h the lower boundary of the model of 
iKoskinen et alj d2007l). The structure is somewhat different 
than those o f lShowman etafl (l2CMl2009h . where there is just 
one jet with the peak approximately located at the 10 _1 bar 
level. The peak flow speed is also much greater in those stud- 
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FIG. 4. — A gravity wave with c = 600 m s , propagating in an atmosphere 
with profiles shown in Figure [5] The horizontal perturbation velocity u ( — 
), intrinsic phase speed, c — uq (•••), and the mean flow acceleration dwo/dr 

( ) are shown. The intrinsic phase speed becomes zero at z/H p fa 5 and 

the wave encounters a critical level. In the layers just below the critical level 
the wave saturates and sheds momentum into the mean flow, causing it to 
accelerate, peaking at a rate over 250 m s rotation - . 

ies at 4 or 5 km s" 1 . It is important to note, however, that these 
differences do not change qualitatively the basic points we are 
making in this paper. 

The Brunt-Vaisala frequency profile N(z) is also shown in 
Figure [3] As already discussed, the profile does not vary 
much over the whole domain: the maximum Brunt-Vaisala 
frequency is just 1.2 times the minimum value. Therefore, 
N does not contribute much to the variation of the index of 
diffraction m. The main contributor to the variation of m is the 
variation of the intrinsic phase speed, which is derived from 
the large variation of flow speeds. This should be compared 
to the analogous terrestrial situation, where the range of flow 
speeds is much lower. This allows N to have a larger effect on 
the variation of m on the Earth. 

3.2. Wave -Background Interaction 
3.2.1. Critical Layer Encounter 

Figure H shows a gravity wave encountering a critical level 
in the upper jet in Figure [5] The wave has a horizontal wave- 
length, 2tt I k, of 2500 km and c = 600 m s" 1 . Here, since 
c > mo as the wave approaches the critical layer from below, 
the momentum deposited in the mean flow causes the mean 
flow to accelerate. This acceleration peaks at over 250 m s" 1 
per rotation. This is large enough to double the flow speed at 
this layer in ^2 rotations — a significant effect. The effect is 
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FIG. 5 . — A gravity wave, with c = -40 m s _1 , propagating in an atmosphere 
described by the profiles in Figure|5] The perturbation to the temperature field 
T ( — ) and the heating rate d7o /df (— • — ) are shown. The wave saturates at 
just above z/H p = 10, where the heating peaks at nearly 80 K rotation -1 . The 
peak energy flux for this wave is approximately 1 W m~ 2 . 

large enough to require its inclusion in a ny simulation of the 
atmospheric circulation (ICho et al.ll2008h . 

The waves encountering critical layers are dissipated. 
Therefore, a flow with a range of flow speeds dissipates all 
gravity waves with phase speeds within this range. That is, 
a spectrum of gravity waves is prevented from propagating 
high into the atmosphere. There are other, secondary effects 
at the critical layer that may also affect the mean flow; but, 
they are not modeled here. They will form the basis of future 
work. For example, the deposition of energy into the flow at 
the critical layer may well lead to the generation of new grav- 
ity waves, which then may propagate further, partly mitigating 
the filtering effect. 

3.2.2. Saturation 

Figure shows an example of a gravity wave saturating in 
the atmosphere of Figure [5] In general, it is possible that a 
wave may not encounter a critical level as it propagates up- 
ward. However, such a wave then travels higher into the at- 
mosphere, where it grows large and, if unabated, eventually 
breaks or suffers dissipation at higher altitudes. Although 
both momentum and energy are deposited in this case, we 
focus here on the effects on the temperature field. 

The wave launched in Figure[5]has c = -40 m s" 1 (i.e., west- 
ward). The horizontal wavelength remains at 2500 km, as in 
the critical layer example of Figure |4] The vertical velocity 
perturbation grows with height. Therefore, so do the zonal 
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FIG. 6. — As in Figure [3] but with the bottom of the computational do- 
main extended down to 100 bars. Here, no is extended downward barotropi- 
cally (independent of height) from the 1 bar level; 7o is extended downward 
so that the profile bel ow the 1 bar level is similar to that of Figure 18 in 
IShowman etaT] <2009l) . 

velocity perturbation u and the potential temperature pertur- 
bation 9, as expected from equation (THY In this case, the 
wave saturates near the top of the jet, where condition ( IT9l ) is 
exactly satisfied. Here, condition (|20T i is approximately sat- 
isfied since the zonal perturbation velocity and the intrinsic 
phase speed (c — «o) are both approximately 122 m s" 1 . The 
saturation deposits energy that causes the atmosphere there to 
heat up. The heating is significant, peaking at ~ 75 K per 
rotation — a 5% change in one rotation. In the absence of 
other effects, the ambient temperature can be doubled in about 
20 planetary rotations (or orbits, assuming 1 : 1 spin-orbit syn- 
chronization). 

A wave that has a phase speed greater than the maximum 
flow speed will not encounter a critical layer. Those with 
phase speeds close to, but still above, the maximum flow 
speed will, in general, saturate in the regions just below the 
maximum flow, as the intrinsic flow speed will be small in that 
region. Similarly, waves with phase speeds just less than the 
minimum flow will saturate as well. In this way the filtering 
effect discussed above is extended beyond those waves with 
phase speeds equal to flow speeds. The main effect of these 
filtered waves on the flow will be lower in the atmosphere. In 
the profile given in Figure [3] where the waves emanate from 
the z/Hp = 1 level, this means that the upper layers of the 
lower jet will be slowed by gravity waves whereas the lower 
levels of the upper jet will be accelerated. 

Those waves that do not dissipate will be able to propa- 
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gate into the upper atmosphere depositing their momentum 
and heat there. Here the changes to the flow can be very large. 
For example, the wave shown in Figure [5] causes a decelera- 
tion of up to 6.8 km s" 1 rotation" 1 as it saturates. This clearly 
dominates the flow at this level. 

Moving the location of the wave origin down does not 
change the basic behavior in the qualitative sense. However, 
the amplitudes are much larger, compared with the case when 
the wave originates higher up in altitude. Thus, the possibility 
exists for stronger effects due to gravity wave interaction with 
the background. 

This is illustrated in Figures [6]and[7] which should be com- 
pared with Figures |3] and[5] Here, we have extended the pro- 
files downwards. The wave is still launched at z/H p = 1, but 
this is now deeper in the atmosphere. The wave again has 
a phase speed of -40 m s -1 , and the horizontal wavelength 
remains at 2500 km. This wave also saturates near the top 
of the upper jet, where the energy deposition into the mean 
flow causes the atmosphere there to heat up. Note that the re- 
gion of heating is lower than when the wave originates higher 
up, as in Figure [7] The heating is significant, peaking at 
~ 3000 K rotation" 1 . The ambient temperature can be dou- 
bled in approximately half of a planetary rotation. In a more 
realistic scenario, dissipation — which we have not included in 
our model — will likely reduce the heating rate. 

3.2.3. Refraction 

So far we have been focusing on the vertical transport of 
momentum and energy by gravity waves. However, the waves 
can als o transport momentum and energy horizontally (cf., 
93.2.51 1. Substituting c = u/k into the the index of refraction, 
equation (fTTT i. and rearranging gives the dispersion relation 
for gravity waves. Here, we consider the case where uq = 0. 
We will examine cases where uq ^0 in the sections following 
this one. 

When there is no background mean flow, the dispersion re- 
lation simplifies to 



LO = ±- 



Nk 



[k 2 + m 2 +l/(4H 2 )] 
We can then use the definitions, 

duj 

U ' = dk 
duo 

dm ' 

to obtain the group velocities. They are: 

N[m 2 +l/(4Hl)] 



1/2' 



U g = ± 



[k 2 + m 2 +l/(4H 2 )] 



3/2 



±- 



Nkm 



[k 2 + m 2 +l/(4H 2 )] 



3/2- 



(23) 

(24a) 
(24b) 

(25a) 
(25b) 



Thus, for propagating waves (i.e., waves for which m is real), 
u g 0. Therefore, gravity waves always propagate obliquely 
and cannot strictly propagate vertically when there is no back- 
ground flow. 

From equations (l25T l we see that -d, the angle of propagation 
with respect to the horizontal, is given by 



tantf 



km 



m 2 +l/(4H 2 )' 



(26) 
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FIG . 7 . — A gravity wave, with c = -40 m s -1 , propagating in an atmosphere 
described by the profiles in Figure|6] The perturbation to the temperature field 
T ( — ) and the heating rate d7o jit (— • — ) are shown. The wave saturates at 
just above z/H p = 13, where the heating rate peaks at just under 3000 K per 
rotation. In terms of the pressure level, this location is actually lower than in 
the case shown in Figuref5] and the magnitude of the peak is nearly 50 times 
greater. The peak energy flux for this wave is nearly 200 W m . Thus, 
having a source at a lower height can have a much stronger effect. 



Since H p is nearly constant, with a value just under 500 km, 
tan# varies with l/m for m > 10~ 6 m" 1 . This gives rise to 
refraction. As a wave propagates into a region of higher m, 
it bends to a more horizontal path. Note that, as we are here 
considering a region with no flow, increasing m is essentially 
equivalent to increasing N. As shown in Figure[3] N increases 
in the thermosphere and so the paths of waves in this region 
will bend towards the horizontal even though the flow is small. 

3.2.4. Trapping 

From equation ( fT2l we can see that in regions where m is 
imaginary, the wave is evanescent: its amplitude decays to- 
wards zero and therefore it does not propagate vertically. This 
can occur when N 2 <Q (i.e., when the atmosphere is convec- 
tively unstable). But, it can also occur for non-hydrostatic 
waves when the buoyancy term becomes dominated by the 
non-hydrostatic term. In addition, when the index of re- 
fraction changes between layers, the wave is reflected at the 
boundary. The amount of reflection is given by the magnitude 
of the coefficient of reflection \r\, where 



m\ —m2 
mi +mi 



(27) 
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FIG. 8. — A gravity wave, with c = 10 m s and horizontal wavelength 
2ir/k — 955 km, trapped in an atmosphere with the structure presented in 
Figure [5] The vertical velocity perturbation w ( — ), mean flow uq ( — ), and 

the real part of the index of refraction m ( ) are shown. The wave is 

trapped in relatively quiescent region between z/Hp K 1.5 and z/H„ m 3.5 
and does not propagate vertically. The region of trapping corresponds to the 
region where m is real. The wave is reflected at the boundaries of this region, 
providing a possibility for resonance. 



In equation d27l) . m\ and ni2 are the indices of refraction in two 
adjacent layers. When ni2 is imaginary, total reflection occurs 
and the wave is evanescent in that region and its amplitude 
decays to zero. However a region of propagation can exist 
between two regions of evanescence. This readily occurs for 
jets, where the intrinsic phase speed varies enough to allow 
the hydrostatic term to dominate in some regions and not in 
others. The region can also occur through variations of the 
Brunt-Vaisala frequency. In Figure [8] we see a wave that is 
trapped in the relatively quiescent region between z/H p « 1 
and z/ Hp w4. Outside this region, the value of Re(m) is small, 
or zero. Trapped in this region the wave is able to interact with 
itself and, under appropriate conditions, resonate. 

This is another mechanism via which waves may be fil- 
tered out by the flow; and so the waves do not reach high 
altitudes at which saturation can occur. However, in this case, 
a trapped wave does not directly interact with a low level flow 
that changes its characteristics. Indeed, between the two re- 
flection layers the wave can propagate horizontally — even in 
the absence of any flow, using the refractive mechanism de- 
scribed above. As long as the layers do not allow much leak- 
age, it is possible for a trapped wave to cover large horizon- 
tal distances — transporting momentum and heat zonally (east- 
west direction). 

3.2.5. Ducting 



As well as being trapped in relatively quiescent regions, it is 
possible for waves to be trapped in a jet. As alluded to above, 
it is possible for such a wave to travel within the region of 
trapping, which is known as a duct or a waveguide. In Fig- 
ure [9] a non-hydrostatic wave with c = 700 m s" 1 is trapped 
within the jet (located at ~5 mbar level) in our model atmo- 
sphere. Note the small values of Re(m) outside the jet. 

In this case, the flows are significant. Therefore, we use the 
full dispersion relation, 



ku + 2H p k [u' + H p (2 (k 2 + m 2 ) u + u' ' ) ] ±2 



OJ = 



H 2 p k 2 



' 1 + 4H 2 (k 2 + m 2 ))N 2 +(u' + H p u[{) 



1/2 



\+4H 2 p {k 2 + i 



(28) 



to develop expressions for the group velocities. 

However, the expressions obtained are large and rather unil- 
luminating. They can be simplified by assuming that u' Q and 
u'q are small. This is realistic since the shear is of the order 
of 10~ 3 s" 1 and u' ' of the order 10~ 8 m" 1 s . This is small 
compared with the other terms in the expressions. This then 
gives 

u g = u + l - p -\^ (29a) 

[k 2 + m 2 +l/(4H 2 )] ' 



[k 2 + m 2 +l/{4H 2 )] yi 

From these expressions, we can see that u g follows uq a s this 
is the larger term on the right hand side of equation (|29al > 
in our model atmosphere. In Figure [9] the values of u g and 



w g are shown. Note that in the center of the duct w g is very 
small while u g is large, so that energy is transported along the 
flow. At the the top and bottom of the duct the vertical group 
velocity increases, while the horizontal group velocity falls. 
Therefore, propagation here is nearly vertical. In Figure [9] 
we show w g as positive, however this is only for the upward 
propagation of energy, at the top of the duct the wave is re- 
flected and the vertical component becomes negative. This 
keeps the wave within the jet. The ray path followed by the 
wave, before the reflection, is shown in Figure [TOl 

The wave can travel large distances in this duct; but, eventu- 
ally, the wave will either escape the duct or be dissipated. The 
range of speeds in the jet may change so that a critical level 
for the wave is created. The wave will then be reabsorbed into 
the flow. Alternatively, if the flow or Brunt-Vaisala frequency 
outside the jet changes so that the buoyancy term is no longer 
small and is dominated by the non-hydrostatic term, then the 
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FIG. 9. — As in Figure[8] but with c = 700 m s and horizontal wavelength 
2ir/k = 1410 km. The horizontal group speed u g ( — ), vertical group speed 

w g (■ ■ ■ ) and the real part of the index of refraction m ( ) are shown. The 

wave, not shown, is trapped in the upper jet between z/H p Ri 3 and z/H p 
9, the region where m is real, and does not propagate vertically above this 
region; it is however, able to propagate along the jet as the large value of u g 
within the jet shows. The wave is reflected at the boundaries of this region, 
providing a possibility for resonance. 

reflection is no longer total and the wave will then leak out 
of the duct. This can, for example, happen when propagating 
into a colder region, assuming the lapse rate remains constant. 
As can be seen from equation (|2"TT i a fall in temperature with 
constant lapse rate will cause an increase in the Brunt-Vaisala 
frequency and thus an increase in the buoyancy term. This 
may occur very far from the original region of wave excita- 
tion. Indeed, in the example given, it is possible to envisage 
jets ducting waves and so transporting energy from the day- 
side of a tidally locked planet to the colder nightside where 
the waves escape the jet and propagate away from the duct 
before dissipating. 

4. IMPLICATIONS FOR CIRCULATION MODELS 

The effects of gravity waves discussed in this work on 
the larger-scale circulation must be parameterized in global 
models because the spatial resolution — both horizontal and 
vertical — required to model them is currently prohibitive. The 
waves important to the large-scale extrasolar planet atmo- 
spheric circulation have horizontal length scales ranging from 
approximately ~ 10 s m to ~ 10 7 m and vertical wavelengths 
as small as ~ 10 4 m. Waves with periods of few hours can 
carry significant momentum and energy flux vertically, but the 
sources of these waves include processes that are not included 
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FIG. 10. — This shows the path of propagation of the wave in Figure[9]as- 
suming that properties of the duct do not change in the je-direction. The path 
shown is the first crossing of the duct, that is until the wave encounters the 
top reflection layer. At this point the wave will reflect and then propagate 
downwards in a mirror image of this path. Note that the wave travels nearly 
one planetary radius before reflection. This means with just six reflections the 
wave will have nearly circumnavigated the planet. Of course in the real situa- 
tion the properties of the duct will change in the ^-direction and the wave will 
probably either leak out of the duct or dissipate before the circumnavigation 
is complete. 

or resolvable in current circulation models. 

The difficulty with representing gravity waves in GCMs ex- 
ists even for the GCMs of the Earth. For example, the param- 
eterization for co nvection is not aime d at producing realis- 
tic gravity waves dCollins et alj|2004h . However, not repre- 
senting gravity waves can affect the accuracy of the GCMs. 
The lack of gravity wave drag can lead to the overestima- 
tion of wi nd speeds, resulting in faster and narrower jets than 
observed dMcLandresslll998l) . Further, gravity waves intro- 
duce UrrbuIejic£_W2flT_™bsecment mixing and thermal trans- 
port dFritts & Alexan der 2003). This leads to greater homog- 
enization of the atmosphere with a reduction in, for example, 
temperature gradients. Gravity waves also interact with plan- 
etary waves, playing a role in important transient phenomena 
(such as sudden stratospheric warming). In the absence of 
gravity waves, thes e phenomena are not accurately modeled 
dRichter et aflEOlOb . 

There are many parametrization schemes currently in- 
corporated or prop osed for general circulation modeling 
dMcLandresslll998l) . In all of the schemes, the basic com- 
ponents are 1) specification of the characteristic of the waves 
at the source level, 2) wave propagation and evolution as a 
function of altitude, and 3) effects on and by the atmosphere. 
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All of them are essentially linear and one-dimensional, in that 
waves only propagate vertically and that only vertical vari- 
ation in the background influence the propagation. As seen 
in this work, linear theory still requires information such as 
the wave's phase speed c and wavenumber£, for example. A 
more complete theory would need spatial and temporal spec- 
tral information. Intermittency is another crucial feature that 
would need to be taken into account. The primary differences 
in various schemes pertain to the treatment of nonlinearity and 
specificity of wave dissipation mechanisms. 

Currently, all global circulation models of hot Jupiters sug- 
gest the presence of a low number of zonal jets. However, 
all the models do not have the resolution required to ade- 
quately resolve gravity waves and are subject to all of the limi- 
tat ions described ab ove. This issue has been previously raised 
by iCho et akl (120081) . in which they advocate caution against 
quantitative interpretation of current model results. For ex- 
ample, without the inclusion of the wave effects discussed in 
this work, high jet speeds and precise e astward shift of puta- 
tive "hot spots" can be questioned (e.g..lCooper & Showmanl 
l200HlKnutson et al.ll2007tlLangton & Laughlinll2007l) ~ 

5. CONCLUSION 

Gravity wave propagation and momentum and energy de- 
position are complicated by the environment in which the 
wave propagates. For example, spatial variability of the back- 
ground winds causes the wave to be refracted, reflected, fo- 
cused, and ducted. Additionally, temporal variability of the 
background winds cause the wave to alter its phase speed. 
Still further complications arise due to the wave's ability to 
generate turbulence, which can modify the source or serve as 



a secondary source, and the wave's interaction with the vorti- 
cal (rotational) mode of the atmosphere. Many of these issues 
are as yet not well-understood and are currently areas of active 
research. 

In this work, we have emphasized only some of these is- 
sues. We have shown that gravity waves propagate and trans- 
port momentum and heat in the atmospheres of hot extrasolar 
planets and that the waves play an important role in the at- 
mosphere. They modify the circulation through exerting ac- 
celerations on the mean flow whenever the wave encounters a 
critical level or saturates. They also transport heat to the upper 
stratosphere and thermosphere, causing significant heating in 
these regions. Moreover, through ducting, they also provide 
a mechanism for transporting heat from the dayside of tidally 
synchronized planets. 

Before relying on GCMs for quantitative descriptions of hot 
extrasolar planet atmospheric circulations, further work needs 
to be performed to ensure that the effects of important sub- 
scale phenomena, such as the gravity waves discussed here, 
are accurately parameterized and included in the GCMs. 
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